Limits to the muon flux from WIMP annihilation in the center of the Earth with the 

AMANDA detector 
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A search for nearly vertical up-going muon-neutrinos from neutralino annihilations in the center 
of the Earth has been performed with the AMANDA-BIO neutrino detector. The data collected 
•i— ( , in 130.1 days of live-time in 1997, ~10 9 events, have been analyzed for this search. No excess 

over the expected atmospheric neutrino background has been observed. An upper limit at 90% 
confidence level has been obtained on the annihilation rate of neutralinos in the center of the Earth, 
as well as the corresponding muon flux limit, both as a function of the neutralino mass in the range 
100 GeV-5000 GeV. 

PACS numbers: 95.35.+d, 95.30.Cq, 11.30.Pb 



I. INTRODUCTION 



There are strong observational indications for the ex- 
istence of dark matter in the universe. Measurements 
of the energy density of the universe, tt , from the com- 
bined analysis of cosmic microwave background radia- 
tion data and high red-shift Type la supernovae favor 
Q = 1, with a matter, fi M , and a cosmological con- 
stant, f^A) component. Combined with data from rota- 
tion curves of galaxies and cluster mass measurements, 
the matter contribution to f2 is 0.3 < fi M < 0.4. Big 
Bang nucleosynthesis calculations of primordial helium, 



lithium and deuterium production, supported by abun- 
dance measurements of these elements, set an upper limit 
on the amount of baryonic matter that can exist in the 
universe, S! B < 0.05 (see Ref. |l| for a recent review of 
values of Q). Non-baryonic dark matter must therefore 
constitute a substantial fraction of Q M . 

In this paper we present results of a search for non- 
baryonic dark matter in the form of weakly interacting 
massive particles (WIMP) using the AMANDA high- 
energy neutrino detector. The next section contains a 
brief motivation for WIMPs as dark matter candidates. 
Section III describes briefly the characteristics of the 



2 



AMANDA detector in the configuration used for this 
analysis. Sections IV and [v| contain a description of the 
simulation and analysis techniques used. In section [vj we 
discuss the sources of the cur rent systematic uncertain- 
ties of our analysis. In section VII we present the results 
of the analysis and we introduce a novel way of calculat- 
ing upper limits in the presence of systematic uncertain- 
ties. An upper limit on the neutrino-induced muon flux 
expected from WIMP annihilation in the center of the 
Earth is derived with this method. A comparison with 
published muon-flux limits obtaine d by e xisting neutrino 
experiments is presented in section VIII. 



II. WIMPS AS DARK MATTER CANDIDATES 

Particle physics provides an interesting dark mat- 
ter candidate as a Weakly Interacting Massive Parti- 
cle (WIMP). The relic density of particle type i de- 
pends on its annihilation cross section, <j, as f2j/i 2 ~ 
3 x 10~ 27 / (crv) (neglecting mass-dependent logarithmic 
corrections), where () indicates thermal average and v 
is the relative velocity of the particles involved in the 
collision (see for example Ref. |J) . Weak interactions pro- 
vide the right annihilation cross section for the WIMPs 
to decouple in the early universe and give a relic density 
within the required range to contribute substantially to 
the energy density of the universe today. This is basically 
what would be needed to solve the dark matter problem. 

In particular, and starting from a completely different 
rationale, the Minimal Supersymmetric extension to the 
Standard Model of particle physics (MSSM) provides a 
promising WIMP candidate in the neutralino, x- The 
neutralino is a linear combination of the B-ino, B, and 
the W-ino, W, the supersymmetric partners of the elec- 
troweak gauge bosons, and of the and Hi], the neutral 
Higgs bosons, and it is stable (assuming R-parity con- 
servation, which is further supported to avoid too rapid 
proton decay) . The actual composition of the neutralino 
can have cosmological consequences since its annihilation 
cross section depends on it. For example, it has been 
argued that a mainly W-ino type neutralino would not 
be cosmologically relevant in the present epoch since it 
would have annihilated too fast in the early universe to 
leave any relevant relic densitycl. 

Still, the large parameter space of minimal supersym- 
metry can be exploited to build realistic models which 
provide relic neutralino densities within the cosmologi- 
cally interesting region of 0.025< f2 x /i 2 <1. Negative re- 
sults from searches for supersymmetry at the LEP accel- 
erator at CERN have set a lower limit on the neutralino 
mass m x > 31 GeV (Ref. ||), while theoretical arguments 
based on the requirement of unitarity set an upper limit 
of 340 TeV (Ref. |||). Imposing in addition the condi- 
tion on fl x h 2 mentioned above, only models with m x < 
10 TeV (Ref. ||) become cosmologically interesting. 

Neutralinos have a non-negligible probability of scat- 
tering off nuclei of ordinary matter. Assuming the dark 



matter in the Galactic halo is (at least partially) com- 
posed of relic neutralinos, elastic interactions of these 
particles with nuclei in the Earth can lead to energy 
losses that bring the neutralino below-,jthe escape veloc- 
ity, becoming gravitationally trappedffij. For high neu- 
tralino masses (> a few hundred GeV) direct capture 
from the halo population by the Earth is kinematically 
suppressed!^. In this case neutralinos can be accreted 
from the population already captured by the solar sys- 
tem. Gravitational capture is expected to result in an ac- 
cumulation of neutralinos around the core of the Earth, 
where they will annihilate. An equilibrium density is 
reached when the capture rate equals the annihilation 
rate. Neutrinos are produced in the decays of the result- 
ing particles, with an energy spectrum extending over 
a wide range of values and bounded from above by the 
neutralino mass. Annihilation of neutralinos directly into 
neutrinos (or light fcrmion pairs in general) is suppressed 
by a factor m^/m 2 due to helicity constraints, where 
is the fermion mass. 

Neutrino detectors can therefore be used to constrain 
the parameter space of supersymmetry by setting liraijts 
on the flux of neutrinos from the center of the EarthnEl 
Note that this indirect neutralino detection will be fa- 
vored for high neutralino masses, since the cross section 
of the resulting neutrinos with ordinary matter scales 
with E„. 



III. THE AMANDA-BIO DETECTOR 

The AMANDA-BIO detector consists of an array of 
302 optical modules deployed in ten vertical strings at 
depths between 1500 m and 2000 m in the South Pole ice 
cap. The strings are arranged in two concentric circles 
of 60 m and 120 m diameter respectively. The modules 
on the four inner strings are separated by 20 m in the 
vertical direction, while in the outer six strings the ver- 
tical separation between modules is 10 m. An optical 
module consists of a photomultiplier tube housed in a 
spherical glass pressure vessel. Coaxial cables (in the in- 
ner 4 strings) and twisted quad cables (in the outer 6 
strings) provide the high voltage to the photomultiplier 
tubes and transmit the signals to the data acquisition 
electronics at the surface. 

Muons from charged-current high-energy neutrino in- 
teractions near the array are detected by the Cherenkov 
light they produce when traversing the ice. The rela- 
tive timing of the Cherenkov photons reaching the opti- 
cal modules allows the reconstruction of the muon track. 
A more detailed description of the detector is given in 
Ref. |l^. The detector was triggered when a majority re- 
quirement was satisfied: an event was recorded if at least 
16 modules had a signal within a predefined time window 
of 2/is. The data taking rate was 100 Hz. 

AMANDA-B10 was in operation during the 1997 
Antarctic winter. The separation of 300 atmospheric 
neutrinos from the data sample collected in that period 
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established the detector as a high-energy neutrino tele- 
scopeEil. The array was upgraded with 122 more modules 
during the antarctic summer 1997/1998 and in 1999/2000 
253 additional ones were added, completing the proposed 
design of 677 optical modules in 19 strings, AMANDA- 

iM 



IV. SIGNAL AND BACKGROUND 
SIMULATIONS 

A. Simulation of neutralino annihilations 

Neutralinos can annihilate pair- wise to, e.g., £ + £~ , qq, 
W+W-, Z°Z°, H? )2 H^, Z°H? >2 and W ± H = F. Neutrinos 
are produced in the decays of these annihilation products. 
Neutrinos produced in quark jets (from e.g. bb or Higgs 
bosons) typically have lower energy than those produced 
from decays of r leptons and gauge bosons. We will refer 
to the first type of annihilation channels as "soft" and to 
the second as "hard" . 

The simulations of the expected neutralino signal were 
done in the framework of the SUSY models described in 
Ref |l3]. The hadronization and decay of therAiinihilation 
products have been simulated using PYTHIAEJ. The sim- 
ulations were performed for six different WIMP masses 
between 10 GeV and 5000 GeV. For each mass, six dif- 
ferent annihilation channels (cc, bb, tt, t + t~, W + W~ 
and Z°Z°) were considered, with 1.25 x 10 6 events gener- 
ated for each. Note that the decay of b- and c-hadrons 
will take place in matter instead of vacuum. This was 
incorporated in the simulations in an effective manner 
justified by the fact that, for the neutralino masses con- 
sidered, the re-interactions of these heavy hadrons with 
the surrounding medium are not dominant, and can be 
parametrized as an effective energy loss at the time of 
decay. As a reference soft spectrum, we chose the anni- 
hilation into bb, and as a reference hard spectrum, the 
annihilation into W + W~. For a given mass, these two 
spectra can be regarded as extreme cases. We have used 
these channels in the analysis described below, bearing 
in mind that a typical spectrum would lie somewhere in 
between. 



B. Simulation of the atmospheric neutrino flux 

Neutrinos from the decay of secondaries produced in 
cosmic ray interactions in the atmosphere constitute the 
physical background to the neutralino search. We have 
simulated this atmospheric neutrino flux using the cal- 
culations of Lipariiij. To obtain the rate of neutrino in- 
teractions producing muons we have used the neutrino 
ancLanti-neutrino— nucleon cross sections from Gandhi et 
alS3. The actual neutrino- nucleon interactions have been 
simulated with PYTHIA using the CTEQ30 parametriza- 
tion of the nucleon structure functions. The use of 
PYTHIA allows to model the hadronic shower produced 



at the vertex of the interaction and, therefore, to calcu- 
late the Cherenkov light produced by secondaries. When 
the neutrino-nucleon interaction occurs within the instru- 
mented volume of the detector, this is a non- negligible 
contribution to the total event light output. 

A three-year equivalent atmospheric neutrino sample 
with energies between 10 GeV and 10 TeV and zenith 
angles between 90° (horizontal) and 180° (vertically 
up-going) has been simulatedEJ. The sample contains 
3.7xl0 7 events, of which 41234 triggered the detector. 



C. Simulation of the atmospheric muon flux 

The majority of the triggers in AMANDA are induced 
by muons produced in cosmic ray interactions in the at- 
mosphere and reaching the detector depth. The simula- 
tion of thiSj-atmospheric muon flux was performed using 
the BASIEV113 program. We note that this program only 
uses protons as primaries. However, the systematic un- 
certainty introduced by this approximation is negligible 
in comparison with that from the present uncertainty 
in the primary flux intensity. Moreover, heavier nuclear 
primaries produce more np*ons per interaction, but with 
lower energies on averaged, which will in general loose 
all their energy and decay before reaching the detector. 
A study performed using the CORSIKAcil air shower gen- 
erator, with the QGSJET option to model the hadronic 
interactions, including the complete cosmic ray composi- 
tion confirms this scenario. 

The simulation of a statistically significant sample 
of atmospheric muon background is an extremely high 
CPU-time consuming task due to the strong rejection 
factors needed. We have simulated 6.3 x 10 10 primary 
interactions, distributed isotropically with zenith angles, 
8, between and 85 degrees, and with energies, E, be- 
tween 1.3 TeV and 1000 TeV, assuming a differential en- 
ergy distribution cx E~ 2 7 (Ref. p2|). The total number of 
triggers produced were 5 x 10 b . Normalizing to the pri- 
mary cosmic ray rate, the generated sample corresponds 
to about 0.6 days of equivalent detector live-time. Due to 
the narrow vertical angular cones used for this analysis 
this background sample is sufficient to model the detec- 
tor response and develop the rejection cuts. In addition, 
a larger sample of background data was used in the train- 
ing of the discriminant analysis program used as cut level 
4. This is described in more detail in the next section. 



D. Muon propagation 

The muons produced in the signal and background 
simulations described above were propagated from the 
production point to the detector taking into account en- 
ergy losses by bremsstrahlung, pair production, photo- 
nuclear interactions and S-ray production from Ref. |2^. 
The Cherenkov light emitted by the secondaries produced 
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FIG. 1: Angular distributions of data and atmospheric muon 
simulation Monte Carlo (MC) at different analysis levels. Top 
to bottom: trigger to level 4. The distributions are normalized 
to the simulated sample, 5xl0 6 events. 



FIG. 2: Angular distribution of muons from atmospheric neu- 
trinos and from the annihilation of neutralinos after filter level 
2. The two extreme neutralino masses and annihilation chan- 
nels considered in this paper are shown. The relative normal- 
ization is arbitrary. 



in these processes is taken into account when calculating 
the response of the detector to the passage of the muon. 



V. DATA ANALYSIS 

The analysis presented in this paper was performed 
on data taken with the 10-string AMANDA detector be- 
tween March and November 1997. The experimental 
data set consists of 1.05xl0 9 events in a total of 130.1 
days of detector live-time. The data were first cleaned of 
noise hits and hits from optical modules that were un- 
stable during the running period. Short pulses, that are 
likely induced by cross talk between channels are also 
rejected at this stage. Details on the data cleaning pro- 

The data is then recon- 



cedure are given in Ref. 24 



structed and five filters consisting of cuts based on the 
event hit pattern and the quality of the reconstruction 
are applied in order to identify potential up-going neu- 
trino candidates. The distributions of the reconstructed 
zenith angle from trigger level (after hit cleaning) until 
filter level 4 for data and simulated atmospheric muons 
are shown in figure [l]. The curves have been normalized 
to the simulated sample, 5xl0 6 events. The uppermost 
curves in the plot show the reconstructed direction with- 
out any quality criteria applied to the fits, showing good 
agreement between the data and the Monte Carlo sam- 
ple along the whole angular range. The curves clearly 
indicate that a small percentage (about 2%) of the origi- 
nally down-going tracks are misreconstructed as up-going 



(cos O less than zero the figure) . The series of cuts de- 
scribed below were developed to reject such misrecon- 
structions, and their effect on the angular distribution is 
also shown in figure for comparison. The filter level 
2 and level 3 curves show that the filtering procedure is 
more effective rejecting the simulated muon background 
than the data. This is due to detector effects not included 
in the simulation of the detector response and surviving 
to these levels, like electronic cross talk between channels 
or inefficiencies of the digitizing electronics. Other pro- 
cesses not included in the background simulations that 
can contribute to the discrepancy are overlapping events 
from uncorrelated cosmic ray interactions and the con- 
tribution from electron neutrino induced cascades. To 
account for this different behavior between data and sim- 
ulated background under standard cuts, we have used an 
iterative discriminant analysis as cut level 4 (see subsec- 
tion VD) trained on a sub-sample of data (which rep- 
resents the real remaining background better than the 
simulations) and a sub-sample of the neutralino signal. 
A final series of high quality cuts were applied after the 
discriminant analysis, bringing the remaining data sam- 
ple to agree with the number of events expected from 
the known atmospheric neutrino flux, as shown in fig- 
ure H and table B. Note that the atmospheric neutrino 
curve and the data curve in figure |]join and follow each 
other in the last two steps of the cuts applied within the 
level 5 filter. The next subsections give a more detailed 
description of the variables used and the cuts applied at 
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FIG. 3: Efficiencies relative to trigger level at filter levels 1 
and 2 as a function of the neutralino mass. 



each filter level. 



A. Filter level 1 

In a first stage, a simple and computationally fast filter 
based on fitting a line to the time pattern of the events, 
was applied to the data sample in order to reject obvi- 
ous down-going tracks. This "line fit" (LF) assumes that 
the known space point of each hit optical module, r^, is 
related to the measured hit time, U, by ?i — r Q + vU. 
The minimization of \ 2 = ~ r a — vU) 2 , where the 

index runs over all the hits in the event, leads to an ex- 
plicit solution for v. The zenith angle of the fitted track 
is readily obtained as cos0 LF = —v z /\v\. The angular 
resolution of the line fit is relatively low since it does 
not incorporate any information about the geometry of 
the Cherenkov cone or about scattering of the Chercnkov 
photons in the ice. Still, its simplicity and computational 
speed makes it a very useful tool for a first assessment 
of the track direction and for rejection of down-going at- 
mospheric niuonsiSI. The first level filter rejected obvious 
down-going atmospheric muons by requiring LF > 50°. 



B. Filter level 2 

The events that pass the level 1 filter are reconstructed 
usiag a maximum likelihood approach (ML) as described 
irJij. In short, the ML technique uses an iterative process 
to maximize the product of the probabilities that the op- 
tical modules receive a signal at the measured times, with 



FIG. 4: Rejection and efficiency at each filter level for the 
data and simulations of the neutralino signal, atmospheric 
neutrinos and atmospheric muons. The dashed part corre- 
sponds to rejection levels surpassing the statistical precision 
of the simulated sample, yielding zero remaining events. The 
neutralino signal curve should be read only with respect to 
the right axis scale, and it shows the relative signal efficiency 
with respect to trigger level. The rest of the curves are plotted 
with respect to the left axis scale. 



the track direction (zenith and azimuth angles) as free 
parameters. The expected time probability distributions 
include the scattering and absorption characteristics of 
the ice as well as the distance and relative orientation of 
the optical module with respect to the tracked. 

The level 2 filter consists of two cuts: the ML- 
reconstructed zenith angle must be larger than 80° and 
at least three hits must be "direct". A hit is defined as 
direct if the time residual, t leB (the difference between the 
measured time and the expected time assuming the pho- 
ton was emitted from the reconstructed track and did 
not suffer any scattering), is small. Unscattered pho- 
tons preserve the timing information. Therefore, the re- 
construction of the direction of tracks with several di- 
rect hits presents a significantly better angular resolu- 
tion. The number of direct hits associated with a track is 
the first quality requirement [applied to the reconstructed 
data and simulated samplesc-i A residual time interval 
between -10 ns and 25 ns was used to classify a hit as 
direct at this level. 

Figure |^ shows the zenith angle distributions of sim- 
ulated muon tracks from neutrinos produced in annihi- 
lation of neutralinos for the two extreme masses used 
in this analysis as compared to that from atmospheric 
neutrinos after filter level 2. The corresponding curve 
for data and simulated atmospheric muons is included in 
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figure y. The combined effect of these two filters on the 
data is a rejection of 98%, as shown in table |. The effi- 
ciencies with respect to trigger level of both level 1 and 
level 2 filters for simulated neutralino signal are shown 
in figure ^, for different neutralino masses and the two 
extreme annihilation channels used. 

Filters 1 and 2 are applied in an initial data reduction 
common to the different subsequent analyses of the data. 
The rest of the cuts described below were specifically 
designed for the WIMP search with the aim of identifying 
and rejecting misreconstructions while maximizing signal 
detection efficiency and background rejections. 

C. Filter level 3 

The angular distribution of the events is the most ob- 
vious difference between the predicted neutralino signal 
and both the atmospheric neutrino flux and the atmo- 
spheric muon background. Neutrinos from neutralino 
annihilations in the center of the Earth would be con- 
centrated in a narrow cone close to the vertical direction, 
while atmospheric neutrinos are distributed isotropically. 
The level 3 filter further restricted the ML-reconstructed 
zenith angle to be larger than 140°, placed a cut on the 
total number of hit modules in the event, N ch > 1 0, and 
on the summed hit probability of the modules with a sig- 
nal, P hit >0.23. The number of hits with time residuals 
between -f ns and 25 ns was required to be larger than 
4 and the number of hits with residuals between -15 ns 
and 75 ns to be larger than 5. At this stage the possible 
correlations between the variables are ignored, and the 
cuts applied to each of them individually. Table Q shows 
the efficiency and rejection power at this cut level. Only 
5xl0 -4 of the simulated atmospheric muon background 
survive this level, compared with 68% of the simulated 
neutralino signal and 10% of the atmospheric neutrinos. 

D. Filter level 4: Iterative Discriminant Analysis 

To account for possible correlations between the 
variables and perform a multidimensional cut in 
parameter space, the next filter level was based 
on an iterative non- linear discriminant analysis, us- 
ing the IDA programed. Given a set of n vari- 
ables, the program builds the "event vector" x fe = 
(xi, ...,x n ,xl,XiX2, ...,xix n ,xl,x 2 X3, where Xi is 

the value of variable i in event k. A class of events, the 
signal or background sample, is characterized by their 
mean vector (x s ) or (xf,), and the mean difference be- 
tween the samples is given by the vector Afi = (x s } — (xf>) . 
The spread of the variables is contained in the vari- 
ance vectors, (j.* = x fc — (x s ) and — x fe — (x;,), 
which are used to define a variance matrix for each class, 

vs ' b = J2k cvts ^ k s ,b(M k ^b) T > where N evts is the num- 
ber of events in the signal or background samples and T 
denotes the transpose. The problem of separating sig- 



nal from background is transformed into the problem of 
finding a hyperplane in event vector space which gives 
minimum local variance for each class and maximum sep- 
aration between classes. This is translated into the re- 
quirement that the ratio R=(a T A^) 2 /a T Va should be 
maximal, where here the variance matrix V is the sum of 
the variance matrices for signal and background and a is 
a vector of coefficients to be determined by training the 
program on a signal and a background sample. A target 
signal efficiency and background rejection factor are cho- 
sen beforehand. The coefficients a are determined in an 
iterative process carried out until the specified rejection 
factor is achieved or a predefined number of iterations 
reached. The coefficients found in this way are used to 
select events from the signal region in the multidimen- 
sional parameter space: each event is characterized by 
the scalar D=a T x and a cut on D serves as the selection 
criterion. 

Eight variables were used in the training of the dis- 
criminant analysis program and in the subsequent cuts: 
the velocity of the line fit, the number of direct hits, the 
number of modules hit, the number of modules hit in the 
string with the largest number of hits, the number of de- 
tector layers with a hit 1 , the extension of the event along 
the three coordinate axes, the average hit probability and 
the probability that the event time pattern is compatible 
with that expected from a vertical up-going muon. This 
set of variables includes combined information from the 
fit track parameters as well as the general spatial and 
temporal topology of the event. 

Since to a first approximation the data consist of at- 
mospheric muon background, seven days of data, evenly 
distributed along the year, were used as the background 
training sample. For the signal training sample, muons 
from the simulations of 250 GcV ncutralinos annihilating 
into a hard spectrum were used. The combination of a 
relatively low neutralino mass and annihilation into the 
hard channel was chosen as giving a "typical" muon spec- 
trum. The target signal efficiency was set to 98% per iter- 
ation and the target global background rejection to 1000. 
The stopping criterion was set to 9 iterations, based on 
the fact that further loops would reduce the number of 
events in the training sample to a too low number to be 
representative of the whole data set. The rejection of 
background achieved was 220 with respect to cut level 3 
since the nine loops were exhausted before reaching the 
desired rejection. The overall signal efficiency attainable 
after the training process is then (0.98) 9 =0.83. The effect 
of the discriminant analysis event selection is shown in 
table |. It indeed achieves the expected signal efficiency, 
retaining 82% of the signal with respect to the previ- 
ous cut level. The discrepancy of the expected number 
of atmospheric neutrinos and the number of remaining 
data events at this level indicates that the data sam- 



The detector was divided in eight horizontal layers of 65 m depth. 
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TABLE I: Rejection of data, of the simulated atmospheric neutrinos and atmospheric-muon background samples and efficiency 
for the simulated neutralino signal from trigger level to filter level 5. 



Filter Level 


Data 
130.1 days 
(events) 


Atmospheric neutrinos 
130.1 days equivalent 
(events) 


Atmospheric muons 
0.6 day equivalent 
(events) 


xx -> ww 

m x =250 GeV 
(% of trigger level) 





1.05xl0 u 


4899 


5x10° 


100 


1+2 


2.3 x 10 7 


2606 


7xl0 4 


79 


3 


1.2 xlO 6 


472 


2588 


68 


4 


5441 


89 


13 


56 


5 


14 


16.0 
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pie is still contaminated by poorly reconstructed down- 
going muons. A last cut level was therefore developed to 
improve the rejection of the remaining misreconstructed 
events and select the truly up-going tracks. 



E. Filter level 5: Final event selection 

The remaining events after the discriminant analysis 
with a zenith angle larger than 165° were passed through 
the following series of cuts. The length spanned by the 
direct hits when projected along the track direction was 
required to be at least 110 m, and the vertical length 
containing all hits was required to be at least 170 m. The 
z component of the center of gravity of the direct hits 
(Zc.o.g. = J2i z i/^dircct hits, where the sum is over all the 
direct hits in the event) was required to be deeper than 
1590 m, and the percentage of hits in the lower half of 
the detector less than 55%. These cuts reject events with 
a spatially uneven concentration of hits, typically due to 
down-going atmospheric muons that pass just outside the 
detector or stop close to the array. 

The remaining data at this level are consistent with 
the expected atmospheric neutrino flux. Figure shows 
the angular distribution of the remaining 14 data events 
and the remaining 16.0 simulated atmospheric neutrino 
events. The angular range shown is for >165°, the re- 
gion where a possible neutralino signal is expected to be 
concentrated. No statistically significant discrepancies 
are found between the expected number of events and 
angular distributions of the atmospheric neutrino back- 
ground and the data. This result is also consistent with 
the results on atmospheric neutrinos presented in Ref. [ll]. 

Due to the different angular shapes of the neutralino 
signal for different neutralino masses (see figure |^ for the 
two extreme cases considered) , we have chosen to restrict 
further in angle the signal region we use to extract the 
limit on an excess muon flux. We use angular cones that 
contain 90% of the signal for a given neutralino mass. 
The remaining data and simulated atmospheric neutrino 
background events for the different angular cones used 
are shown in table 0. The background rejection power 
and signal efficiency from filter level 1 to 5 are shown in 
figure |] along with the effect on the data sample. 



VI. SYSTEMATIC UNCERTAINTIES 

An essential quantity when deriving limits, as we do 
in the next section, is the effective volume, V at{ , of the 
detector. It is the measure of the efficiency to a given 
signal and it is defined as 



where n L5 is the number of signal events after filter level 
5 and n gon the number of events simulated in a volume 
V sm surrounding the detector. The effective volume of 
AMANDA-B10 as a function of muon energy is shown 
in figure [?]. Given a MSSM model producing a muon 
flux with a given muon energy spectrum, the effective 
volume of the detector for this particular signal is also 
calculated through equation [l|. This is shown in figure |§| 
for the different neutralino masses used in this analysis. 
The shaded bands in both figures indicate the systematic 
uncertainty estimated as described below. 

The evaluation of V Btc is subject to experimental and 
theoretical systematic uncertainties present in the anal- 
ysis. We have performed a detailed study of the effect of 
the uncertainty in several variables on the resulting ef- 
fective volume by propagating variations in any of them 
to the final evaluation of V bS . 

Measurements of the scattering and absorption 
lengths, A s and A a , using pulsed and DC light sources 
deployed with the detector at different depths and YAG 
laser light sent from the surface through optical fibers, 
have shown that these quantities exhibit a depth de- 
pendence which is correlated with dust concentration at 
different levels in the iceEE A simulation of the detec- 
tor response including layers of ice with different optical 
properties has been developed and used to evaluate its 
effect on the results. The effects introduced are muon- 
energy dependent and therefore dependent on the neu- 
tralino model. The effective volumes calculated with the 
layered ice model are reduced between 1% and 20% with 
respect to the uniform ice model, except for the lower 
neutralino mass and soft annihilation channel (100 GeV) 
where the effect reaches 50%. 

A further correction accounts for the uncertainties in 
the optical modules' total and angular sensitivities. It 
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FIG. 5: Angular distribution of the remaining data events 
(dots) and simulated atmospheric neutrino events (shaded 
area) at filter level 5. The angular range shown is between 
f65° and 180°. The shaded area represents the total un- 
certainty in the expected number of events. 
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FIG. 6: Angular distribution of the remaining fraction of 
neutralinos at filter level5 with respect to trigger level from 
the two extreme neutralino masses studied in this paper. 
The angular range shown is between 165° and f80°. 



is known that during the process of re-freezing after de- 
ployment, air bubbles appear in the column of ice that 
has been melted, changing locally the scattering length 
of the ice and distorting the effective optical module an- 
gular sensitivity with respect to that measured in the 
laboratory. We have used a specific ice model for the 
ice in the holes that accommodates this effect. The fact 
that it appears after deployment and that it is not di- 
rectly measurable in the laboratory makes it difficult to 
assess. Only by an iterative process of comparison of 
data and different hole-ice models can it be quantified. 
We estimate this effect to yield and increase of 20% in 
effective volume with respect to the uniform angular re- 
sponse model with, again, the soft annihilation channel 
of the lowest mass neutralino giving a stronger effect of 
34%. An additional 20% uncertainty on the total optical 
module sensitivity has been used. 

The way to combine of all these effects into a final esti- 
mate of the total uncertainty in V ett is a difficult subject, 
since they are not independent contributions. As de- 
scribed in the previous paragraphs, by varying the initial 
parameters used in the simulations of the detector and 
in the ice properties, we have obtained a range of possi- 
ble values for the effective volume, which we consider as 
equally probable giving our current understanding of the 
detector. We have chosen to take the nominal V elf to be 
used in equation [l] as the middle value of this range. As 
a conservative estimate of the uncertainty we take half 
the width of the range of values obtained. We thus con- 
clude that our current estimate of V et { is affected by a 



systematic uncertainty crv eS /V ea between 10% and 25% 
depending on the neutralino mass considered, the lower 
mass of 100 GeV giving the larger relative error. A sim- 
ilar estimate including the same effects has been made 
for the atmospheric neutrino Monte Carlo. In this case 
we estimate the uncertainty on the effective volume for 
atmospheric neutrinos to be 20%. 



Further uncertainty in the number of expected atmo- 
spheric neutrinos (column 3 in table |) is caused by 
the uncertainties present in the calculation of the atmo- 
spheric neutrino flux. This is estimated to be of the order 
of 30% in the energy region relevant to this analysis, and 
originates mainly from uncertainties in the normalization 
of the primary cosmic pay spectrum and in the hadronic 
cross sections involved^. This has been taken into ac- 
count as an additional effect on top of the experimental 
uncertainty on the effective vo lume f or atmospheric neu- 
trinos, as described in section VII B. 



It has recently been shown that different muon prop- 
agation codes can produce differences in the muon flux 
and energy spectrum at the detector depth (see for ex- 
ample Re^pil). The code used in this analysis uses the 
LohmannEj parametrizations for muon energy loss, which 
produce results in agreement within about 10% of more 
recent codesE3 for muon energies up to a few of TeV. 
We have not included any systematics arising from the 
treatment of muon propagation in the ice in this analysis. 



9 



TABLE II: Number of data events, simulated atmospheric 
neutrino background events and the corresponding N 90 for 
the angular cones containing 90% of the signal for the dif- 
ferent neutralino masses. These angular cuts are applied in 
addition to the level 5 filter described in section |^. "s" and 
"h" denote the soft and hard annihilation channels. The num- 
bers in parenthesis in column 5 show N go obtained without 
including systematic uncertainties. 



niv 


Angular cut 


Data 


Atmospheric 


Non 
1 1 90 


(GeV) 


(deg) 


(evts.) 


neutrinos (evts.) 




100s 


167.5 


10 


12.1 


9.2(4.7) 


lOOh 


168.5 


9 


10.8 


6.6(4.7) 


250s 


170.0 


7 


8.6 


5.9(4.1) 


250h 1 
500s / 


172.0 


5 


6.1 


5.6(3.9) 


1000s 


173.0 


4 


4.6 


5.3(3.9) 


500h 


173.5 


4 


4.6 


5.3(3.9) 


lOOOh 1 








5.6(4.7) 


3000s J 


174.0 


4 


3.9 


3000h s 










5000s \ 


174.5 


3 


3.9 


4.4(3.6) 


5000h > 











VII. RESULTS 

From the observed number of events, n obs , and the 
number of expected atmospheric neutrino background 
events, n B , an upper limit on the signal, N^, at a cho- 
sen confidence level /3%, can be obtained. We have used, 
the unified approach for confidence belt construction^ 
to calculate 90% confidence level limits. In section B be- 
low we briefly describe a novel way of calculating limits 
in the presence of systematic uncertainties that we have 
used to obtain the final numbers presented in this paper. 



A. Flux limits: the standard approach 

For detectors with a fixed geometrical area A, it is 
natural to derive a muon flux limit directly through <f>^ < 
N^/A ■ t, where t is the detector live-time. However, due 
to the large volume of AMANDA and the lack of sharp 
geometrical boundaries it is the effective volume V eti , as 
defined in equation |l|, that has to be used to determine 
a limit on the volumetric neutrino-to-muon conversion 
rate, T vll . The effective volume provides a measure of 
the detector efficiency since, in addition to through-going 
tracks, it takes into account the effect of tracks starting 
or stopping within the detector. A limit can then be set 
on r„ M , that is, on the number of muons with an energy 
above the detector threshold E tl „ produced by neutrino 
interactions per unit volume and time, 



N 90 



(2) 



K« • t 

T vtl includes all the detector threshold effects and 



model dependencies, as indicated below, and can be di- 
rectly related to a more physically meaningful quantity, 



the annihilation rate, 
the Earth through 



r„ M (m x ) = 



r A , of neutralinos in the center of 



4tt#| jo 



2^ B xx^x 



dE v 



(3) 



x a 



(E V \E^ > E thI ) p N dE v , 



where the term inside the integral takes into account the 
production of muons through the neutrino-nucleon cross 
section, ov+at, weighted by the different branching ratios 
of the XX annihilation process and the corresponding neu- 
trino energy spectra, B x ^x dN v jdE v . p N is the nucleon 
density of the ice and is the radius of the Earth. We 
have used a muon energy threshold of 10 GeV in the sim- 
ulations of the signal, which has been taken into account 
through the muon production cross section. 

Equation [?] is solved for T A . T A depends on the MSSM 
model assumptions, as well as the galactic halo model 
used, being related to the capture rate of neutralinos in 
the Earth. Different neutralino models predict different 
capture and annihilation rates that can be probed by ex- 
perimental limits set on T A . The right column of table III 
shows the limits thus derived for T A . The correspond- 
ing curves are shown in figure ^. Quoting limits on the 
annihilation rate has the advantage that the detector ef- 
ficiency and threshold are included through equation |^ 
and, therefore, numbers published by different experi- 
ments are directly comparable. This is not usually the 
case when presenting limits on muon fluxes, where at 
least the detector energy threshold enters in a non triv- 
ial way and prevents direct comparison between experi- 
ments. However, since it is common in the literature to 
present limits on the muon flux per unit area and time, 
we transform below our limit on T A into a limit on the 
muon flux from neutralino annihilations in the center of 
the Earth. 

The total number of muons per unit area and time 
above any energy threshold E thr within a cone of half 
angle 9 C as a function of the annihilation rate is 



>E thr ,e>e c ) = 



dE„ 



d e- d2N » 



dEf.de ' 



(4) 



where the term d 2 N fl /dE fl d0 represents the number of 
muons per unit angle and energy produced from the neu- 
tralino annihilations, and includes all the MSSM model 
dependencies for neutrino production from neutralino an- 
nihilation and the neutrino-nucleon interaction kinemat- 
ics, as well as muon energy losses from the production 
point to the detector. The upper limits on the annihila- 
tion rate are thus converted to a limit on the neutralino- 
induced muon flux at any depth and above any chosen 
energy threshold and angular aperture. The 90% confi- 
dence level upper limits on the annihilation rate and the 
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TABLE III: The 90% confidence level upper limits on the muon flux from neutralino annihilations in the center of the Earth, 
<ji> M , for a muon energy threshold > 1 GeV. The last column shows the threshold-independent neutralino annihilation rate, F A . 
Detector systematic uncertainties have been included in the calculation of the limits. The corresponding numbers without 
including uncertainties are shown in parenthesis 



m x [GeV] 


annihil. 




r A 




channel 


[xl0 3 Km- 2 yr _1 ] 




100 


hard 


8.9 (6.3) 


4.0 (2.9) xlO 14 




soft 


133.5 (68.2) 


4.3 (2.2) xlO 16 


250 


hard 


2.1 (1.5) 


1.3 (0.9) x 10" 




soft 


6.9 (3.9) 


3.8 (2.2) xlO 14 


500 


hard 


1.5 (1.1) 


2.5 (1.8) XlO 1 * 




soft 


2.7 (1.9) 


4.4 (3.0) xlO 13 


1000 


hard 


1.5 (1.2) 


6.5 (5.4) xlO 11 




soft 


1.8 (1.4) 


9.2 (6.8) xlO 12 


3000 


hard 


1.1 (1.0) 


7.5 (6.7)xlO iu 




soft 


1.5 (1.3) 


1.5 (1.3)xl0 12 


5000 


hard 


1.1 (1.0) 


3.2 (2.8)xlO iu 




soft 


1.5 (1.2) 


7.6 (6.4) xlO 11 



muon flux at an energy threshold of 1 GeV derived using 
equations |2|, || and [I] are shown in parenthesis in table [I] . 
The fluxes have been corrected for the inefficiency intro- 
duced by using angular cones that include 90% of the sig- 
nal, so the numbers presented represent the limit on the 
total muon flux for each neutralino model. The thresh- 
old of 1 GeV has been chosen to be able to compare with 
published limits by other ex perim ents that have similar 
muon thresholds (see section VIII| below). 



B. Evaluation of the limits including systematic 
uncertainties 



However, the best limits an experiment can set are af- 
fected by the systematic uncertainties entering the anal- 
ysis. Including the known theoretical and experimental 
systematic uncertainties in the calculation of a flux limit 
is not straightforward, and often overlooked in the lit- 
erature. A precise evaluation of a limit should involve 
the incorporation of both the uncertainties in the back- 
ground counts, cr b , and in the effective volume, a v . An 
additional caveat arises since the uncertainty in the effec- 
tive volume introduces in turn an additional uncertainty 
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FIG. 9: 90% confidence level upper limits on the neutralino 
annihilation rate, r A , in the center of the Earth as a func- 
tion of the neutralino mass and for the two extreme anni- 
hilation channels considered in the analysis. The dashed 
lines indicate the limits obtained without including sys- 
tematic uncertainties and correspond to the numbers in 
parentheses in table The symbols indicate the masses 
used in the analysis. Lines are to guide the eye. 



FIG. 10: 90% confidence level upper limits on the muon 
flux at the surface of the Earth, MI as a function of the neu- 
tralino mass and for the two extreme annihilation channels 
considered in the analysis. The dashed lines indicate the 
limits obtained without including systematic uncertainties 
and correspond to the numbers in parentheses in table [II 
The symbols indicate the masses used in the analysis 
Lines are to guide the eye. 



in the expected number of background events, on top 
of the 30% uncertainty used in the background neutrino 
flux <j b . A proper implementation of the systematics in 
the calculation of a limit should take this correlation into 
account. 

One approach to incorporate systematic uncertainties 
into an upper limit has been proposed in Ref. |34|. We 
have developed a similar method suited to our specific 
case which includes the systematic uncertainty in V c{{ in 
the calculation of N 90 used in equation |^. The method-is 
a modified Neyman-type confidence belt constructional 
The confidence belt for a desired confidence level j3 is 
constructed in the usual way by integrating the Poisson 
distribution with mean n tot =n s +n B so as to include a (3% 
probability content. But the number of events for signal 
and background, n s and n B , are taken themselves to be 
random variables obtained from Gaussian distributions 
with means equal to the actual number of signal and 
background events observed and widths corresponding to 
the systematic uncertainties in signal and background. 

Given an experimentally observed number of events, 
N oxp , the 90% confidence level limit on the number of 
signal events is obtained by simply inverting the calcu- 
lated N 90 (n tot ) at the corresponding n tot =N oxp value. In 
this way the different uncertainties for signal and back- 
ground and the correlation between them are included 
naturally. 



In summary, the inclusion of our present systematic 
uncertainties in the flux limit calculation yields results 
which are weakened between ~10% and ~40% (practi- 
cally a factor of 2 for the soft channel of m x = 100 GeV) 
with respect to those obtained using Ngo calculated with- 
out systematics. The effect is dependent on the WIMP 
mass, and it reflects the better sensitivity of AMANDA 
for higher neutrino energies. Figures | and show the 
90% confidence level limit on the neutralino annihilation 
rate and the corresponding limit on the resulting muon 
flux for a muon threshold of 1 GeV for the hard and soft 
annihilation channels considered in the analysis. The 
symbols show the particular neutralino masses used in 
the simulation. The lines are to guide the eye and they 
show the limits obtained including systematic uncertain- 
ties (solid). The dashed lines, included for comparison, 
show the values obtained using the Neyman construc- 
tion with the unified ordering scheme without including 
uncertainties. Table III summarizes the corresponding 
numbers. 



C. Effect of neutrino oscillations 

To account for neutrino oscillations among the differ- 
ent flavors, the atmospheric neutrino spectrum should 
be weighted by a factor W(E 1/ ) which includes the prob- 
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FIG. 11: The AMANDA limits on the muon flux from neu- 
tralino annihilations from figure [H^ compared with published 
limits from MACRO, Baksan and Super-Kamiokande. The 
dots represent model predictions from the MSSM, calculated 
with the DarkSUSY packageH. The dashed area shows the 
modelsrdiisfavored by direct searches from the DAMA collab- 
orations. 



ability that a muon neutrino has oscillated into an- 
other flavor in its way through the Earth to the detec- 
tor. For the purpose of illustration consider a two-flavor 
oscillation scenario between and v T . Then W(E„) 
= 1 - sin 2 (2(9) sin 2 (l.27 Am 2 [eV 2 ] D^kmJ/E^GeV]), 
where is the diameter of the Earth, 9 the mixing 
angle and Am 2 the difference of the squares of the flavor 
masses. Note that the effect depends strongly on the neu- 
trino energy and it is negligible in the high energy tail of 
the atmospheric spectrum since the oscillation length is 
then much larger than the Earth diameter. If we choose 
sin 2 (2#)=l and Am 2 = 2.5 x 10" 3 eV 2 based on the re- 
sults obtained in Ref. |3^, the number of expected atmo- 
spheric neutrino events is reduced between 5% and 10%, 
depending on the angular cone considered. This would 
weaken the limits by about the same amount. 

The effect of neutrino oscillations on the possible 
WIMP signal is model dependent and has been estimated 
in Refs. ^ and |3^. However the authors reach different 
conclusions on the direction of the effect: up to a fac- 
tor of two in increased muon flux in Ref. 37 and a re- 
duction of about 25% in Ref. B$ for a neutralino mass 



of 100 GeV. For higher neutralino masses both authors 
predict a less pronounced effect, which, becomes negligible 
for the higher masses considered ir£3 (m x > 300 GeV). 
We have not included any oscillation effect on the neu- 
trinos from the WIMP signals considered in this paper. 



VIII. COMPARISON WITH OTHER 
EXPERIMENTS AND THEORETICAL MODELS 

Searches for a neutrino signal from WIMP annihila- 
tion in the center of the Earth have been performed by 
MACRO, Baikal, Baksan, and Super-J^amiokande-. 

In figure |l| the jupsults of BaksanH, MACRO0 and 
Super-KamiokandcEj are shown along with the limits 
from AMANDA obtained in the previous section and 
theoretical predictions of the MSSM as a function of 
WIMP mass. In order to be able to compare with the 
other experiments, the Super-Kamiokande limits have 
been scaled by a factor 1/0.9 to represent total flux lim- 
its, instead of limits based on angular cones including 
90% of the signal as originally presented in Ref. [l2| The 
90% confidence level muon flux limits for a muon energy 
threshold of 10 GeV published by the Baikal collabora- 
tion range between 0.63xl0 4 km~ 2 yr _1 for a zenith half 
cone of 15° and 0.54xl0 4 km -2 yr _1 for a zenith half 
cone of 5° (Ref. |]|). Since these results are not pre- 
sented as a function of WIMP mass, and are quoted at 
a slightly higher muon energy threshold, we have not in- 
cluded them in the figure but we mention them here for 
completeness. 

Each point in the figure represents a flux obtained 
with a particular combination of MSSM parameters, fol- 
lowing Ref. ^J. The sixty four original free parame- 
ters of the general MSSM have been reduced to seven 
by the standard assumptions about the behavior of the 
theory at the GUT scale and about the supersymme- 
try breaking parameters in the s-fermion sector. The 
independent parameters left are the Higgsino mass pa- 
rameter /i, the ratio of the Higgs vacuum expectation 
values tan/3, the gaugino mass parameter M 2 , the mass 
m A of the CP-odd Higgs boson and the quantities m D , A t 
and A b from the ansatz on the scale of supersymmetry 
breaking. These parameters were varied in the follow- 
ing ranges: -5000< fi <5000 GeV , -5000< M 2 <5000 
GeV, 1.2< tan/3 < 50, m A < 1000 GeV, 100< m D < 3000 
GeV, -3m Q < A b < 3m Q and -3m Q < A t < 3m D . Mod- 
els based on parameters already excluded by accelerator 
limits are not shown, and the figure is restricted to those 
models which give cosmologically interesting neutralino 
relic densities, 0.025< fl x h 2 <0.5. A local dark matter 
density of 0.3 GeV/cm 3 has been assumed. Theoretical 
predictions for high mass neutralino models lie below the 
scale of the plot, since in this case the number density of 
neutralinos falls down rapidly if the dark matter density 
is kept fixed. 

A complementary way to search for neutralinos is by 
measuring the nuclear recoil in elastic neujtralino-nucleus 
collisions on an adequate target material]. Experiments 
using this direct detection technique set limits on the 
neutralino-nucleon cross section as a function of neu- 
tralino mass. The same scan over MSSM parameter space 
used to generate the theoretical points in figure [ll] can 
be used to identify parameter combinations that are ac- 
cessible by direct searches. There is not, however, a one- 
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to-one correspondence between the results of the direct 
detection searches and the expected neutrino flux from 
the models probed, so comparisons with the results of in- 
direct searches have to be performed with care. We have 
indicated the models disfavored by the DAMA collabo- 
rations by the dashed area in the figure, which has to be 
taken as an approximate region in view of the mentioned 
difficulties in comparing both types of detection tech- 
niques. We note that the models that yield high muon 
fluxes, and that are disfavored by both current results 
from direct searches and by the limits shown in the fig- 
ure, have in common a low value of the mass, around 
92 GeV. 



IX. SUMMARY 

We have performed a search for a statistically sig- 
nificant excess of vertically up-going muons with the 
AMANDA neutrino detector as a signature for neutralino 
annihilation in the center of the Earth. Limits on the 
neutralino annihilation rate have been derived from the 
non-observation of a signal excess over the predicted at- 
mospheric neutrino background. We have included the 
effect of the detector systematic uncertainties and the 
theoretical uncertainty in the expected number of back- 
ground events in the derivation of the limits, presenting 
in this way realistic limit values. 



A comparison with the results of MACRO, Super-K 
and Baksan, as well as with theoretical expectations from 
the MSSM are presented. AMANDA, with only 130.1 
days of effective exposure in 1997, has reached a sensitiv- 
ity in the high neutralino mass (>500 GeV) region com- 
parable to that achieved by detectors with much longer 
live-times. 
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